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Abstract 

Heat transfer results have been obtained for randomly varying laminar flow 
in a heated parallel plate channel. Two types of heat transfer boundary con- 
ditions are considered: (a) The channel walls are at constant temperature; 

or (b) a specified uniform heat flux is transferred at the walls. The heat 
transfer behavior is obtained along the entire length of the channel. The 
temperature distributions are two-dimensional but a slug flow velocity pro- 
file is assumed. The results show significantly higher mean wall tempera- 
tures and, for short channels, lower wall heat fluxes when compared to the 
steady flow cases. 


INTRODUCTION 


Although flows in internal flow systems are commonly assumed to be steady, 
velocity fluctuations are always present to some extent. These flow fluctua- 
tions may arise due to flow instabilities or from unsteadiness in the pumping 
system. It is of interest to determine the effect of the random velocity fluc- 
tuations on such things as wall heat transfer and wall temperatures which will 
also fluctuate and may be significantly different from the steady flow case. 

Practically no analytical papers have dealt with the problem of heat transfer 
in a parallel plate channel with randomly fluctuating velocities. In Ref. 1 
an analysis was carried out for a sinusoidally varying laminar slug flow in 
a channel. For the constant wall temperature boundary condition and constant 
properties the time average wall heat transfer was found not to change appre- 
ciably due to the oscillating flow velocity. These results were extended to 
the case for the two-dimensional velocity profile in Ref. 2 and similar re- 
sults were found. 

In Ref. 3 the heat conduction through a solid with randomly varying internal 
heat generation or randomly varying wall heat fluxes and wall temperature was 
analyzed. Since the heat conduction problem is linear a straightforward 
analysis yielded analytical solutions. 

The present analysis treats the problem of heat transfer to a fluid with ran- 
domly varying velocity. Since this problem is nonlinear, it is difficult to 
find analytical solutions as in the linear case. The present method of 
analysis depends partially on model sampling or the Monte Carlo technique to 
obtain the solution to the present problem. The randomly varying velocity- 
time history is generated by randomly choosing velocity amplitudes and times 
between changes in velocity amplitude from the appropriate distributions. 

The velocity function so generated can be used in obtaining the wall heat 
transfer or temperature as a function of time using standard analytical 
methods. By repeating this process so as to obtain ensemble averages the 
various characteristics of interest of the system can be obtained. 
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Random Velocity 

The fluid velocity is assumed constant across the channel but varying in 
time. Nondimens lonall zing the velocity by the mean flow we can write 

U(9) - 1 + yV'(6) (l) 


The U* is the fluctuating component of velocity with a mean value of zero. 
The value of U* is constant for a time increment Q§ and then U' changes 
to a new value. The amplitude factor V is a parameter of the problem and 
determines the magnitude of the velocity fluctuation. 


The times between changes of the velocity amplitude (6g) are assumed inde- 
pendent events randomly distributed along the time axis with an average time 
of occurrence 6 & . The Poisson distribution (Refs. 4 and 5) 


r 6 / e t 


f(e,K) 


(e/ejK 


kT 


(2) 


gives the distribution of K random events occurring independently during 
time 6 for a given average time per event of 0 & . The occurrence of one 
event at the end of time 0g is given by first no event occurring during 
time ©g and then one event occurring during the incremental time 68. The 
probability of this occurrence is given by the product of the two probabili- 
ties 

f(e 8 ) - f(e B ,o)f(ae,i) - e“ e ^ e& (3) 

a 

The distribution of 6§, the time between velocity changes, is thus given by 
f(©g,) . To pick a time between change in velocity from the distribution 
f(0g) we proceed as shown in Ref. 4. Set R, the random number picked from 
a uniform distribution between zero and 1, equal to the cumulative distribu- 
tion of f (8). 

“%/^a 

f(©)d© = 1 - e (4) 

Thus by having the computer generate R using a uniform random number gen- 
erator we can solve for the random variable 0g, the time between changes in 
velocity from Eq. (4a). We can rewrite the above, using R * 1 - R, as 

G t “ " e a ln R Ub) 

After finding 0g the time to the change in value of U* a new value of 
U* must be randomly chosen from f(U* ), the distribution of amplitudes of 
the fluctuating velocity component U* . In the present analysis we assume 
that U* takes the values of +1 or -1 with equal probability. The new 
value of U* is thus chosen at the end of time increment 0g by generating 
a new value of R. If R is less than 0.5, then U' is taken at +1 
otherwise it is -1. In this manner an entire history of U can be gener- 
ated (see Fig. 2). Using a similar procedure different distributions of 
amplitude of U' can be used as well as different distribution of time to 
change 0g to represent different randomly fluctuating velocities. 

The value of X taken in the present analysis is always less than 1 so 
that the fluid velocity U is always positive. Since the average value of 
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\J' is zero, (U' ) ■ 1 X (l/2) - 1 X Cl/2) - 0, then the mean value of U is 
given by (U) ■ 1 + r(U') - 1. The atandard deviation of the velocity is 
given by 

a„ - (<u2) - <u) 2 ) l/2 . r (5) 

Since 

(u 2 ) » <(i + ru') 2 > - i + r^u'S) + 2r<u*> - 1 + r 2 


The final term of interest is the normalized autocovariance of the velocity 


0, ,(A) = (WMe**)) - we)) 2 

<u 2 (s)> - <u(e )> 2 


( 6 ) 


where A is a time Increment added to 6 . 


If we let f(0,O) be the probability that an event has not occurred in time 
0 in Eq. (2), then since there are only two possibilities, either an event 
has not occurred or has occurred in time 0, we can write 


(u(e)u(e + a)) = (u 2 )f ( a,o) + <u> 2 [i - f(A,o)l 


Then 


Cjj(A), the normalized autocovariance, becomes 


Cu(e) 


-e/e 

e a ((u 2 > - (u> 2 ) -e/e 

■ 1 ■ 1 o e & 


(?) 


( 8 ) 


It can be seen that in the present case the mean, (U) = 1, the mean deviation, 

-e/e_ 

cry » r, and the normalized autocovariance, Cy * e , contain all the pa- 
rameters needed to describe the fluctuating velocity (Eqs. (l) and (3)). 
Knowing these quantities allow the generation of randomly fluctuating veloci- 
ties by model sampling that have the correct statistical behavior as was done 
in the present case. 


Analysis for Uniform Wall Temperature 

The unsteady energy equation for flow in a parallel plate channel is (Fig. l) 
dt , f s dt d 2 t / Q * 

57 + u(t) “ a rr (9a) 

dy 


Viscous dissipation and axial conduction have been neglected and radial con- 
vection is assumed zero. Constant properties are assumed. The slug flow 
approximation has been made for the velocity. This equation for the uniform 
wall temperature case can be written in dimensionless form as 



+ u(e) 


dT _ d^T 

3x £y2 


(9b) 


vising as the dimensionless temperature T = (t - t Q )/(t v - t Q ) (see Fig. l) ; 
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dimensionless axial distance from entrance of heated channel X - (xx)/(Uj_a 2 ) 
and dimensionless time 6 - (ta/a 2 ); where represents the mean flow in 
the channel. 

The boundary conditions on Eq. (9) are 

T ■ 0 at X ■ 0 for all 6 and Y, constant entering temperature (10a) 

T ■ 1 at Y ■ ±1 for all X and 9, constant wall temperature (10b) 

^ ■ 0 at Y = 0 for all a and 9, symmetry (10c) 

The method of solution foDAows the procedure given in Ref. 1. The solution 
for steady laminar slug flow in a channel with a constant wall temperature 
is given by 


T_ 


1-2 




cos I^Y 


n=0 


( 11 ) 


where E- are the eigenvalues [n + (l/2)]n. 


For the case of randomly fluctuating velocity we assume a solution, as in 
Ref. 1, of the form 



( 12 ) 


This form of the solution automatically satisfies the boundary conditions 
(10b) and (10c). Substitution of Eq. (12) into Eq. (9) gives 

dG n _p , * 

35* + u ST " <°n < 13 > 


Using the method of characteristics we write the auxiliary equations as 




(14) 


The first two terms in Eq. (14) can be integrated to yield a set of curves 
or characteristic lines in the x - 9 plane. 



(15) 


The characteristic lines of interest begin at X * 0 since it is at this 
point the boundary condition (10a) must be satisfied. Equating the first 
and last terms of Eq. (14) yields a relationship for G n as a function 
of 9 along the characteristic line. 
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( 17 ) 


The boundary condition 0 n * 1 at 0 - 0 fixes the value of CL at 
X *f 0 since 0 - 0 Q corresponds to the beginning of a characteristic line 
at X ■ 0. This fulfills boundary condition (10a). Substitution of Eq. (17) 
into Eq. (12) gives 


T(0) - 1 - 2 


^2 igi 

n«0 


(18) 


In a situation where the wall temperature is specified it is desirable to 
compute the heat transferred from the wall to the fluid. This can be found 
from the temperature distribution by applying Fourier's law 


<1 



(19) 


Differentiating Eq. (18) 
flux 


Q(0) 


qa 

k(t- t } 


and evaluating it at 

-4< e - e o> 

n**0 



Y 


1 yields the wall heat 


( 20 ) 


To solve for Q(0) in Eq. (20) 
can generate a random velocity 
can be used to give values of 


✓*0 

x(e) = J u(e)d9 


we have to find 0 - 0 Q from Eq. (15). We 
U(0) as discussed previously. This U(0) 

X as a function of 0 where 


This is plotted in Fig. 3. We would like to obtain values of the wall heat 
flux Q at various values of time, Q(0), Q(0 + A), . . . Q(6 + x&) for a 
fixed position along the heated channel X. For Q(0) the value of (0 - 0 Q ) 
needed in Eq. (20) is obtained as shown in Fig. 3. A value of 0 is picked. 
This defines a point along the X coordinate. Moving back along the X co- 
ordinate a distance X defines point 0~ as shown graphically in Fig. 3. 

To find Q(0 + A) we need the value of [(8 + A) - 0 Q ^] given by the follow- 
ing equation ' 



This is found as before by finding a point along X in Fig. 3 corresponding 
to 0 + A moving back a distance X along X then finding 0 , . In this 
manner we can obtain ^Q(0). ^Q(0 + A), lQ(0 + 2A) . . . based on* the randomly 
fluctuating velocity lu(0) where the 1 superscript denotes one member of 
an ensemble. This same procedure is repeated using a newly generated func- 
tion of U(0) for the second member of the ensemble by generating a new ran- 
dom velocity 2 U(0) and finding 2 Q(0), 2 Q(0 + A), . . .. This is done K 
times and the following ensemble averaged values are obtained; <J(0), the 
mean wall heat flux; the standard deviation of the wall heat flux; • 


Sq - {[tfw] - [d(S )] 2 } l/2 


( 22 ) 
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and the normalized autocovariance of the vail heat flux 

<*(*> . (» 

These results have been evaluated and vlll be given later. The value of K 
was taken as 2000 after comparison of results for a K of 4000 shoved neg- 
ligible change. 


Analysis for Uniform Wall Heat Flux Case 
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For this case the energy equation is the same as given in Eq. (9) except the 
dimensionless temperature is now T* « (t - t Q )k/qa. The boundary conditions 
remain the same but Eq. (10c) becomes 

dH* 


W 


» ±1 at Y « ±1 for all X and 0 


(24) 


Following the analysis in Ref. 1, the steady slug flow solution for constant 
wall heat flux is given by 


Tj = X + - 2 ^ ^ e"^ nn ^ X cos nffY 

For the present analysis we assume a solution of the form 


(25) 


T* = F o (0X) + 


3Y 2 - 




LHP. r (ex) cos nxY 
(nn) 2 “ 


(26) 


Substitution into the energy equation gives 

aF ° * „ dF ° , 
dF_ aF„ , 

sr + u sr ’ - (IW) F n 

Using the method of characteristic* the auxiliary equations are 
dX 


&e = 


u 


dF 


de 


ax 


dF 


(27a) 
(27b) 

(28a) 
(28b) 

(29a) 

F n = e (29b) 

where 9 - 0 is given by Eq. (14). We can now write the solution for the 
wall temperature by letting Y = 1 


u -(njr) 2 F 


These can be solved to give 
1 " 6 o 

-(nrt) 2 (e-e_) 


F o - 0 - e o 


l 
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(6 




1 -(n*) 2 (e-0 o ) 

- — — m 

( n *) 2 


(30) 


Then following the procedure described earlier we can then find ^T£(0), 
!tJ(0 + A), !t*( 0 + 2A), . . . based on the randomly fluctuating velocity 
3-Ul (0). This procedure is repeated and the following ensemble averages are 
obtained. The mean wall tenq?eratur e , T£(0); the standard deviation of the 
wall temperature, 


S T - ([^(a) 2 ] - [^(0)] 2 } 1//2 (31) 

and the normalized autocovariance of the wall temperature 

0 w * ■*> - (&»>]* 


These result? have been evaluated and will be given later. 


Limiting Solutions for Large Times Between Changes in Velocity 


For very large values of 0 a , the average time between changes in velocity, the 
problem reduces to the case of heat transfer in s ,eady flow, half of the time 
at velocity U « 1 + y and half the time at U « 1 - y. This means that the 
wall temperature and wall heat transfer results will be given by the steady 
flow results half the time based on velocity U ■ 1 + y and the other half of 
the time based on U = 1 » y. In Fig. 4 is a plot of the steady wall temper- 
ature versus X for the case of uniform wall heat flux. The wall temperature 
at the steady flow velocity U *» 1 + y can be read off this curve at 


x i+r = 


XX 

(1 + yju^a 2 


X 

1 + r 


(33) 


So that to find the wall temperature at X when the velocity is U » 1 + y 
read the steady flow wall temperature at X^^ which is equal to X divided 
by 1 + y. For the limiting case of & & «* » the average wall temperature for 
the uniform wall heat flux case is now given by 



(* 5 , 1-1 



(34) 


We can find the limiting case for the fully developed region as follows. In 
the fully developed region for the steady wall heat flux case the steady wall 
temperature is given by T* 8 => X + (l/3) . We can then rewrite Eq. (34) for 
the fully developed regions' as 


rrrfc 

W,00,d 




1 

3 


(35) 


We can see from Eq. (35) that in the limiting case of large times between 
velocity changes in the fully developed region that for zero value of. the 
fluctuating velocity, y = 0 the mean wall temperature is the same as the 
steady state case, however, when the velocity fluctuation approaches the 
value of the mean flow, y + 1 so that the velocity will be close to zero 
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port of the time then the mean wall temperature becomes very large and ap- 
proaches infinity. The wall temperature standard deviation for the constant 
wall heat flux case for the limiting case of very large time between velocity 
changes is given following Eq. (31) es 


®T,w " [^l+r ” 2T * ^,1-r 



(36) 


This reduces to 

s Tj „ - ~ ts7) 

We can find the limiting result in the fully developed region which is 



(38) 


We can see that the standard deviation of the wall temperature for 6 very 
large in the fully developed region goes to zero as the amplitude of the ve- 
locity fluctuation T go to zero. However, as Y -* 1 then S T w d -* «. 


A similar type of limiting argument would apply to the wall heat flux Q for 
the constant wall temperature boundary condition. As in Eq. (34) we can write 
for the case of very long times between velocity changes, 6 -* », 



(39) 


where and are evaluated from the steady flow result as shown in 

Fig. 6 for Q at x 1+y . or x-j - as given by Eq. (33), We can similarly 
write the standard deviation of the wall heat flux for tt:e limiting case as 
in Eq. (37) 



Q i+y ~ Qj-r 

2 


These limiting values are discussed in the section on results. 


Results for Uniform Wall Heat Flux Case 


The wall temperature results for the uniform wall heat flux case are shown in 
Fig. 4. The solid line shows the limiting solution of the wall temperature 
for the steady flow case. Mean wall temperatures are shown for two different 
amplitudes of flow velocity fluctuations r “ 0.9 and y » 0.5. The Y of 
0.5 means the fluctuating component of the velocity has a positive or negative 
amplitude 50 percent of the mean flow. These mean wall temperature results 
are considerably above the steady flow values. The larger value of X giving 
results further above the steady flow case. The mean wall temperature results 
are also given for different values of 6 & , the average time between changes 
in the fluctuations in the velocity. These results indicate that for large 
values of 9 & the average value of the mean wall temperature fall further 
above the steady flow result. The limiting result of 0 a -* <» are also 
shown. This limiting solution is the average of the wall temperatures for 
the steady flow case when the velocities are taken at the maximum value 
1 + Y*, and the minimum value 1 - y. 

Also shown in Fig. 4 are the standard deviations of the wall temperature S^. 
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These results indicate the amplitude of the variation of the instantaneous 
wall temperature above and below the average wall temperature. The standard 
deviation * is given in Fig. 4 as the difference in values between the open and 
the solid symbols. The standard deviation approaches zero as the average wall 
temperature approaches the steady flow value. Very large values of the stan- 
dard deviation which indicate large fluctuations in amplitude of the wall tem- 
perature occur for large values of y and 9 & as can be seen in the figure. 

The wall temperature autocovariance results are shown in Fig. 5. For compari- 
son the autocovariance for the velocity fluctuations are shown as a solid line. 
The slope of the autocovariance curves are indicative of the average time be- 
tween changes in the wall temperature, the steeper the slope the shorter the 
time between changes in wall temperature. In Fig. 5 the slopes are smaller 
for the unsteady wall temperature then for the unsteady velocity, meaning that 
the average time between changes in wall temperature are somewhat longer than 
for the case of the fluid velocity. Also for larger values of X the average 
time between changes in wall temperature increase as can be deduced from 
Fig. 5. 

Using the mean value, the standard deviation and the autocovariance results 
for the wall temperature, the wall temperature can be generated as a function 
of time in a similar manner as the velocity used in the present analysis was 
generated as discussed earlier. 



Results for the Uniform Wall Temperature Case 


The resulting mean wall heat flux for the constant wall temperature boundary 
condition is shown in Fig. 6. The solid line shows the resultxng wall heat 
flux for the steady flow case. Mean wall heat fluxes are shown for two differ- 
ent amplitudes of the velocity fluctuations X *» 0.9 and X = 0.5. The mean 
wall heat flux is significantly below the steady wall heat flux for heated 
channel lengths less than about X of 0.6. However, since the fluid must 
finally reach the temperature of the wall in the case of constant wall temper- 
ature boundary condition the mean wall heat flux for values of X greater 
than 0.6 must be above the steady value. These results indicate that for 
short heated channels there can be a significant decrease in heat transfer for 
unsteady flow. The mean wall heat flux results are also shown for different 
values of 9 & , the average time between changes in the velocity fluctuations. 
The results indicate that for larger values of X and larger values of 0 
the mean wall heat flux falls further below the steady flow value for short 
channels. The limiting results for 9 & °° are also shown. Also shown 

in Fig. 6 are the standard deviations of the wall heat flux Sq. These 
results are given as the difference between the solid and the open points. 

The standard deviation becomes smaller as the mean wall heat flux comes 
closer to the steady flow case, i.e., small values of X and 9 & . 


The wall heat flux autocovariance results are shown in Fig. 7. For comparison 
the autocovariance for the velocity is shown as a solid line. The slopes of 
the heat flux results are smaller than for the velocity indicating a larger 
average time between changes in wall heat flux compared to the average time 
between changes in the velocity fluctuations. The slope was smaller for 
larger values of X than for smaller values of X. 
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Conclusion 

• 

Adding an unsteady axial component of velocity with a mean value of zero to 
slug flow In a heated channel can result In significantly higher mean vail 
temperature and, for short channels, lover wall *eat fluxes. These effects 
were found to increase as the amplitude of the velocity fluctuations were in- 
creased and also as the average time between velocity fluctuations were in- 
creased. 

The method of analysis consisted of generating a randomly fluctuating velocity 
as a function of time by picking amplitudes and times to change in amplitudes 
randomly from the appropriate distributions. These were used in the energy 
equation which were solved using standard methods of solution. This process 
was repeated so as to obtain ensemble averages. This method allows complex 
nonlinear stochastic problems to be readily solved that would be very difficult 
using more usual analytical procedures. The resulting mean value, standard 
deviation and autocovariance results for wall temperature and wall heat flux 
given in the present report can be used as in the case of the randomly fluc- 
tuating velocity to generate an approximate time history for the wall temper- 
ature and wall heat flux that could be used in the analysis of other problems 
such as temperature stress in materials. 

The present method of analysis can be extended to different forms of the ran- 
domly fluctuating velocity. 

Several assumptions were made that could have significant effects on the re- 
sults including: no random fluctuations in velocity normal to the mean flow 

direction, no axial conduction and the slug flow velocity profile. Further 
analysis removing these restrictions would be of interest. 
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Figure 2. - Fluid velocity. 
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Figure 4. - Wall temperature mean, T w ; and standard 
deviation, Sj; results for constant wall heat flux case. 
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Figure 5. - Normalized auto covariance of wall temperature for constant 
wall heat flux case, y - 0.9. 
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